Voltage induced conformational changes and current control 
in charge transfer through molecules 
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Transport through molecular contacts with a sluggish intramolecular vibrational mode strongly 
coupled to excess charges is studied far from equilibrium. A Born-Oppenheimer approximation in 
steady state reveals voltage dependent energy surfaces which cause abrupt conformational changes 
of the molecular backbone. These are directly related to transitions between current plateaus which 
are relatively robust against thermal fluctuations. In a regime accessible in experiments this allows 
to operate a molecular junction as a current switch or as a molecular machine in form of a valve 
controlled by time dependent bias and gate voltages. 

PACS numbers: 85.65.+h,73.23.Hk,73.63.-b,85.85.+j 



I. INTRODUCTION 

In recent years impressive progress has been achieved 
in fabricating, functionalizingand controlling molecular 
junctions for charge transfer [lH3| and as nano-machines 
[J]. This includes advanced techniques to link molecu- 
lar structures to metallic electrodes and gates, to design 
molecular rectifiers or to induce rotary molecular motion. 
Particularly fascinating is the interplay of charge transfer 
and conformational dynamics. 

Theoretically, density functional calculations (DFT) in 
combination with non-equilibrium Greens function tech- 
niques as well as model based descriptions in terms of 
master equations have elucidated fundamental aspects of 
single molecule contacts. A particular challenge is to cap- 
ture electromechanical properties of these devices where 
injected charges drive internal molecular degrees of free- 
dom such as vibrational or rotational modes far from 
equilibrium which in turn back-act on the transport chan- 
nels. A well-established procedure is to determine energy 
surfaces for these modes from the conformation of the 
molecular backbone in or at least very close to equilib- 
rium. Since couplings between molecules and electrodes 
are typically weak and relevant intra-molecular excita- 
tions are fast compared to displacements of the skeleton 
of the entire structure, this approach has been success- 
fully applied to explain a variety of experimental observa- 
tions [3|. What happens, however, far from equilibrium 
is much less understood. 

To address this issue we consider in this paper a min- 
imal system consisting of a conjugated organic molecule 
with two sub-units the relative orientation of which de- 
termines its transfer properties when contacted to elec- 
trodes. A typical realization is e.g. a biphenyl compound 
with the dihedral angle between the two benzene rings 
acting as conformational degree of freedom. This set- 
up has received considerable attention in recent exper- 
iments, see e.g. 5 8]. So far, however, only situations 
close to equilibrium (low voltage regime) and for very 
weak electrode-molecule couplings have been addressed. 
Theoretical work has revealed the role of vibrational exci- 



tations and dynamical symmetry breaking [Icj within 
Master equation approaches. For compounds with fixed 
torsion mode DFT calculations apply to show micro- 
scopic details of conductance properties [Tll - [l3j . 

Here we focus on the complementary regime of strong 
non-equilibrium with a rotational mode that is not 
frozen, but may approach a voltage dependent optimal 
configuration in steady state. This scenario requires 
electrode-molecule couplings large compared to vibra- 
tional but small compared to electronic excitations on 
the molecule (details are given below). Then, a Born- 
Oppenheimer approximation (BO) applies which pro- 
vides voltage dependent energy surfaces. Notably, stable 
molecular conformations can be controlled by bias and 
gate voltage which are directly related to a robust quan- 
tization of current transfer. This opens the possibility for 
realizations of molecular machines as e.g. switches 

and valves. 

For charge transfer through molecular contacts BOs 
in steady state have been discussed previously in several 
contexts. In [l5| a type of self-consistent treatment for 
a contact with a single electronic level coupled to single 
harmonic intramolecular degree of freedom revealed the 
appearance of hysteresis and negative differential con- 
ductance. Greens functions methods have been used in 
conjunction with a BO to extract effective potential sur- 
faces for nearly harmonic vibronic modes in [jq . How- 
ever, both treatments exploit that the driven dynamics 
of an harmonic oscillator is analytically known. Instead, 
in the sequel a methodology for substantial anharmonic- 
ity of intramolecular modes and/or Coulomb interaction 
between excess charges on several molecular electronic 
levels is studied. 



II. MODEL 

The biphenyl molecule (Fig. [I} consists of two benzene 
rings which can be twisted by an angle <f>. Each of the 
rings may carry one excess electron and the tunnel cou- 
pling between left ring (L) and right ring (R) depends on 
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FIG. 1: Sketch of the biphenyl molecule showing the rota- 
tional angle <j>. Connection to the electrodes is established 
via the thiole groups. 



leads. For large Stark shift we instead get 9 — and the 
eigenstates become the localized ones that only couple to 
a single lead each. 

In eigenstate representation the molecule Hamiltonian 

is 



Ho = -^^ + £o(0)|O)<O|+i? + (<«l+}<+l 
+E-(<f>)\-)(-\+E D (4,)\D)(D\. 



(4) 



this torsion angle as the dominant intramolecular mode 
with respect to transport_properties. The total molecule 
Hamiltonian then reads 



H 



V(^)+T^)(d[d R + d R d L ) 
Vs, 



+E C (n L + n R ) + U Q n L n R + -^{nL - n R ) ,(1) 



where the first line describes the dynamics of the torsion 
angle and tunneling between the rings, i.e., 

V((f>) = F cos 2 (2^), T{<j>) = T cos</> , (2) 

and the second line captures the charging energy E c 
for individual excess charges on the molecule and the 
Coulomb repulsion Uo for double occupancy, respectively. 
Following microscopic calculations [Tljj we also incorpo- 
rate a Stark energy shift Vs of the molecular sites due 
to an external electric field. It turns out that this Stark 
shift has substantial impact on the molecular conforma- 
tion and transport properties. Further we have the local 
electron densities n a = d^ a d a ,a = L,R, where operators 
d' L (d R ) create excess electrons on the left (right) sub- 
unit. Crucial for the charge transfer across this structure 
is the separation of energy scales with a small height 
Vq = 0.05eV for the rotational barrier of the neutral 
molecule, a large hopping element T = 0.5eV, and a 
large inertial moment I = 20000eV _1 p7lfl8| . As it is 
typical for conjugated organic molecules, orbitals are thus 
delocalized throughout the molecular backbone. Accord- 
ingly, it is appropriate to introduce the molecule's charge 
energy eigenstates as linear combinations of the left and 



right localized states: d} + = cos 9d 



svn9d} L + cos6d R . The mixing angle 9 G [— 7r/2,7r/2] 



sm.9d R and d 



depends on the torsional degree of freedom 



= arctan 



T{4>) 



(3) 



with the eigenenergies of the one-electron Hamiltonian 
\{4>) — ^JT(4>) 2 + Vgji. As one can see, for vanishing 
Stark shift, we have 9 — ±7r/4 (the sign of 9 jumps at 
(j) = 7r where the hopping matrix elemnt crosses zero) and 
the energy eigenstates become the symmetric and anti- 
symmetric combinations of localized states, which couple 
symmetrically (except for a sign) to the right and left 



Here, |0) denotes the neutral molecule, |±) the two one- 
particle eigenstates, and \D) is the doubly occupied state. 
The torsional dependent energies of these states are given 
by: E {cj>) = V(<f>), E±(ct>) = V(<f>) ± \(<f>) + E c , and 
E D (cf>) = V(cf>) + 2E c + U . 

The molecular junction is now modeled by the Hamil- 
tonian H = Hq + Hi + Hl + H R with H a = J2k( e ka ~ 
Ha)c], a Cka,ct £ {L, R} describing left (L) and right (R) 
lead, respectively, as reservoirs of non-interacting quasi- 
particles with creation (annihilation) operators c\ a (cfc Q ). 
The chemical potentials ^l/ R = V g ± Vt/2 are fixed by 
the bias voltage Vb and the gate voltage V g . The Stark- 
shift varies with the bias according to Vs = k Vb with a 
parameter k < 1. Coupling of the left (right) lead to the 
left (right) ring is described by 



Hi = jLd L tp L + l R d R ip R + h.c. 
= 7i(cos^<it_ — sm.6d_)il>L 

+7H (sin 9d} + + cos9d^)4> R + h.c. 



(5) 



with tjj a = ^2 k c a ^k- Charge transfer into/out of the 
molecule is determined by rates T a = D a j^/2Ti where 
electrodes are taken in the wide band limit with a den- 
sity of states D a . Apparently, the coupling of the molec- 
ular states depends on the mixing angle 9 which itself 
depends sensitively on the Stark-field according to ([3]). 
For instance, the coupling of the |+) state to the left 
(right) lead is given by 2r L cos 2 (6») [2r fl sin 2 (0)]. As 
shown in Fig. [21 for Tl — T R = T these couplings are 
identical only for vanishing electric field, i.e. Vs = and 
9(4>) = 7r/4. For finite values of Vs this symmetry is im- 
mediately broken, e.g. around <fi = n/2 — S with small 
deviations S, one has 9(ir/2 — 8) « (To/Vs)5 so that 
2rcos(6>) 2 « 2r > 2rsin 2 6» « 2T(T 5/V S ) 2 . 

In the remainder we are particularly interested in 
the steady state current (I a ) — eTr{WoN a (t — > oo)} = 
eTr{N a W(t — > oo)} where N a = J2 k k c a ^k is the num- 
ber operator in lead a and W(t) is the density operator of 
the junction. With electrodes residing in thermal equilib- 
rium the reduced density operator p{t) = TrL ea ds{l / l / (0} 
provides the relevant transport information. A standard 
procedure is to derive from the Liouville-von Neuman 
equation p(t) = (— i/?i) Triads [-ff, a master equa- 

tion by treating Hi as a perturbation (weak coupling). 
This way, upon employing the usual Born-Markov ap- 
proximation one arrives at p(t) — (—i/h)[Ho, p(t)]+Rp(t) 
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charge transfer through the contact suggests an alterna- 
tive approach though which avoids this deficiency. In the 
spirit of a BO one first calculates steady-state solutions 
/? s t(</>) for the electronic system at any fixed value of the 
rotational angle and then uses this solution to extract an 
effective steady-state potential for the torsional mode. 
Within the basis {|0), |±), \D)} and for given <f> corre- 
sponding energies are not degenerate [apart for Vs = 
at 4> — j mod7r (Fig. [3J]. A secular approximation 
applied for fixed (f> to © leads thus in steady state to 
R{poo,P++,P — ,PDd)U = with 



FIG. 2: Frank-Condon-factor cos 6((f>) as a function of (j> and 
Vs. As described in the text, a finite electric field breaks 
the symmetry of the one-electron eigenstates (6(<p) = 7r/4) 
and induces a preference for angles <f> ~ n/2 in the coupled 
system. 



with the Redfield tensor (HH^ 

Mt) = -i / dTTr Lea ds[#/, [Hj{-t), W(t)]} . (6) 
n Jo 

Positively definite steady-state solutions are only ob- 
tained if a secular approximation is applied to ([5]) so 
that couplings between populations (diagonal elements) 
and coherences (off-diagonal elements) of non-degenerate 
states are dropped [21|. The corresponding Redfield 
tensor R determines the steady state density due to 
Rpst = 0. Further, from© one has 



{la) 



-i lim / dTTr{W(t)[I a ,H!(-T)]} 



(7) 



with the current operator I a — (ie/h)[N a , Hi] = 
{ieja/h) Y^k(^ aC ka ~ ^ aC ko)- *~*f course, in steady state 
In 1 1, 



III. STEADY STATE OBSERVABLES 

For low voltages (|I4|, \ V g \ < Vq) and very weak cou- 
plings, transport properties are restricted to the low en- 
ergy sector. Thus a representation of p(t) in the ba- 
sis {|0, m), | — , n)} is appropriate, with vibrational modes 
{|m), \n)} probing only the vicinity of (j> = on surfaces 
E ((f>),E-((f>), respectively 0113] (see Fig.©. 

In regime of higher voltages, however, this scheme 
poses severe difficulties. On surfaces Eo,E± the density 
of torsional eigenstates (typical level spacings \/Vq/I ~ 
meV) strongly increases such that they should more ac- 
curately be described in terms of coherent states. This 
necessitates the inclusion of coherences in ^ also for 
non- degenerate eigenstates which typically leads to neg- 
ative populations and thus to a breakdown of the pertur- 
bative treatment. A separation of time scales between 
the sluggish motion of the dihedral angle and a faster 
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Here, transition rates read 000 



Vout 
^±0 



and odd = 
self-energies are given by 

/ j n ab,a 



Vout 
^D+ 



(8) 



= Sip + S m , a ±± = 
S^J 1 !, where incoming 



*5> 



,b, a f a [E a (<P) - E b (0)] (9) 



and outgoing ones by £°f = hJ2 a ^ab, a - S^ Q . Fur- 
ther, we have introduced Fermi distributions f a (E) = 
1/{1 + cxp[/3(i? — and coupling rates: r o b iQ = 

2T a \(a\d a + d^|6)| 2 . As noted above, these coupling 
rates are proportional to sin 2 9(4>) and cos 2 #(</>), respec- 
tively, and reduce to r Q in absence of a Stark-shift 
(V s = : 6 = tt/4). 

The angle dependent steady state density determines 
all relevant quantities such as the Born-Oppenheimer sur- 
face (BOS) for the torsional mode 



v eS (<t>) = e MpjM 

je{o,+,-,D} 



(10) 



(red line in Fig. [3]) and the angle dependent current 
through lead a [see (O] 



-poo(£+o :Ce 

-P-^-to 



/V°ut v^ n ~\ 



out 

D + , 



^D-M) 



(red line in Fig. rjj Eventually, the Schrodinger equa- 
tion for the torsional degree of freedom with the effective 
potential (fl"0|) is solved providing eigenfunctions 
and energies e v . Note that these depend on the elec- 
tronic steady state of the junction and thus on bias and 
gate voltage as well as temperature. Angle-averaged 
expectation values of observables X(<f>) are obtained 
from (A(0)> = A (0)| (0)| 2 exp(-/3e„)/Z 

with partition function Z. A thermal equilibrium state 
of the torsional sector is not in conflict with steady-state 
solutions of the electronic sector far from equilibrium due 
to the coupling of the former one to a heat bath environ- 
ment of residual vibronic modes. Note that this in turn 
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FIG. 3: Energy surfaces of the bare molecular states |0), |±) 
and the BOS 14ff vs. the torsional angle <fi for 14 = 0.5V, 
V g = 0.5V . In upper part bare molecular eigenstates are 
used, in the lower one a Stark shift between the left and right 
site of Vs = 0.314 is assumed. 



suppresses any possible bistability of the system, e.g. in 
Fig. El (top) a bistability between the global minimum sit- 
uated at cf> = and the local one around <f> = tt/2. With 
typical energy barriers in V e g of order eV (Fig. [3]), the 
mean angle (</>} is dominated over a broad temperature 
range mostly by torsional states well localized around 
the minima of V e s (green lines in Fig. Finite angle- 
averaged currents (/(</>)) result from sufficient overlap of 
angle dependent currents (fTTj) with these states. 

Before we proceed, let us specify the constraints for 
the above scenario. It is based on a time scale separation 
between electronic passage through the molecule and tor- 
sional motion (adiabaticity) and on a separation between 
electronic transition energies and electrode-molecule cou- 
plings hT (perturbation theory) . Specifically, this means 
e„ <C hT <C E + — E- which should be accessible in ex- 
perimental set-ups due to e„ ~ meV <C E+ — £L ~ eV. 



FIG. 4: Angle dependent current and probability density for 
the T — 80K rotational state. Parameters are the same as in 
Fig. [3] 




FIG. 5: Net current in units of eT (top) and mean torsional 
angle (degrees, bottom) in steady state vs. gate voltage V g 
and bias voltage 14 (in units of V) with Vs = 0. 



IV. TRANSPORT PROPERTIES 

Let us now discuss the steady state properties in more 
detail. We start with the case of a vanishing intramolec- 
ular field k — 0. In Fig. [5] (bottom) the twist angle 
displays stability plateaus which reveal conformational 
changes with increasing bias and gate voltage, respec- 
tively. In a central diamond (low bias voltage, V g > 0) 
the molecule resides in an almost planar configuration 
((4>) I~k w 0), while outside it is found mainly out of plane 
with {(f>) /tt ~ j or (</>) /?r« |. One additional plateau ap- 



pears with (4>) /tt w when four transition channels are 
accessible: |0) — > |±), |±) — > \D). This step-wise snap- 
ping leaves its imprint on the current as seen in Figs. [S] 
(top) and |H According to (fTTj) and (HJ) transport chan- 
nels open when eV exceeds transition energies between 
molecular states which lie in the conduction window de- 
termined by V g . The current exhibits distinct plateaus 
with (I)/(eT) = 0, |, g, |, 1 which can be individually ac- 
tivated by the gate voltage. An asymmetric response to 
V g is also seen with a central Coulomb diamond shifted 
towards positive voltages [loj . 
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FIG. 6: Current and equilibrium angle vs. V g at fixed Vt = 
3y,for k — 0. Arrows denote the opening (closing) of transi- 
tion channels. 




In realistic junctions the presence of an intramolecular 
field (according to |ll| we set k = 0.3) breaks the sym- 
metry, thus opening a gap between the E + and the £L 
surfaces (Fig. EJ) and decreasing the HOMO (£ , )-LUMO 
(E- ) separation [ll[ . While for the mean twist angle and 
the net current (Fig. [7]) the general appearance of the cen- 
tral diamond survives, new wedge-like structures appear 
with (4>)/tt = 1/2 and thus (I) » due to T(tt/2) = 0. 
Namely, for finite Stark shifts the effective coupling be- 
tween the leads and the one-particle eigenstates becomes 
highly asymmetrical and peaked near <f> — n/2. This ef- 
fect is surprisingly robust and dominates the conduction 
properties even at very low internal fields (k = 0.01). 
The fully symmetric situation for Vs = turns thus out 
to be extremely unstable and not applicable to actual 
junctions. A cut through the Vb-V g plane at a fixed 
bias voltage (Fig. HJ) reveals the details of the opening 
and closing of transitions between molecular states upon 
sweeping V g from negative to positive values and the si- 
multaneous snapping of the dihedral angle. In contrast, 
previous treatments where the torsion angle is consid- 
ered to be fixed [6l48Ml2| provide conductances oc cos 2 (0). 
Here, transitions between current plateaus occur rather 
abruptly with simultaneous switchings of the molecular 
backbone. Notably, this behavior is very insensitive to 
temperature fluctuations; the system showed just the 
same plateaus (albeit slightly rounded) for T — 300K. Of 
particular interest are sequences (I)/eT w 1 f> H 1/2 
in the regime V g — ... 4 with substantial conforma- 
tional changes ((/>)/ n wl/4f> 1/2 o- 0. This could be 
exploited as a molecular switch (current on/off) with the 
benefit of being insensitive against variations of tempera- 
ture and internal fields in contrast to alternative propos- 
als. In another application the voltage dependence of the 
molecular conformation could be used by operating the 
biphenyl as a molecular motor [l4j |. More specifically, in 
a first version a slowly varying time periodic gate voltage 
at a fixed bias voltage induces a valve-like behavior (e.g. 
open for cf> = 0, closed for <f> = ir/2). The state of the 
rotor is directly read off from the variations of the dc- 
current. In an extended set-up side-groups are attached 



FIG. 7: Same as Fig. [5] but with V s = 0.314. The vertical 
line at VJ, = 2>V marks the cut displayed in Fig. [8] 
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FIG. 8: Current (I) in units of eF (solid) and average twist 
angle (cf>)/iv (dashed) at V b = 3V vs. V g (in V). To the 
left (right) of the maximal current arrows indicate the open- 
ing (closing) of transition channels between molecular states 
|0),|±>,P). 



to each of the biphenyl rings to break in an additional 
static electric field the 7r-symmetry. A time-dependent 
gate voltage may then generate left- resp. right-handed 
rotations of the individual rings. This situation will be 
explored in a subsequent work. 



V. SUMMARY 

We have analyzed charge transfer through conjugated 
organic molecules contacted to metallic electrodes where 
a single sluggish vibrational degree of freedom strongly 
couples to excess charges. In the regime of higher volt- 
ages and molecule-electrode couplings large compared 
to vibrational excitations, but small compared to in- 
tramolecular electronic level spacings, voltage dependent 
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energy surfaces are determined. Bias and gate voltage is gratefully acknowledged, 
allow to activate current plateaus that correspond to spe- 
cific molecular conformations under steady state condi- 
tions. Applications may include non-linear junctions ro- 
bust against thermal fluctuations and molecular motors. 
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